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ABSTRACT 


Explicit  formulas  are  given  for  the  error  at  an  arbitrary  point 
in  the  electron  density,  spin  density  and  form  factor  derived  from  an 
approximate  wave  function.  In  the  derivation  extensive  use  is  being  made 
of  some  previous  results  by  Kinoshita.  The  main  treatment  is  preceeded 
by  a  short  discussion  of  some  of  the  more  mathematical  aspects  of  the 
problem  and  a  simple  numerical  example  is  given  in  the  last  paragraph. 
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1.  INTRODUCTION 


The  properties  of  an  atom  in  its  ground  state  can  be  calculated  from  the 
solution  of  the  time  independent  SchrOdinger  equation 

HV.  >  £„T<. 

where  Eq  is  the  lowest  eigenvalue  of  the  Hamilton  operator  H  and  4*0  is  the 

corresponding  eigenfunction  which  depends  on  the  space  and  spin  co-ordinates 

of  the  electrons.  Unfortunately,  however,  the  SchrOdinger  equation  is  of  such  a 

complicated  form  that,  except  for  the  simplest  case  of  the  hydrogen  atom,  it 

does  not  seem  possible  to  obtain  the  exact  solution.  On  the  other  hand  there  are 

methods  available  by  which  approximations  to  vpQ  of  varying  degree  of  accuracy 

can  be  obtained.  As  more  and  more  refined  calculations  had  been  made,  at  least 

on  small  atoms,  there  has  been  a  continued  interest  in  deriving  limits  of  error 

1)  2) 

for  the  expectation  value  of  the  energy  '  and  other  quantities  '  calculated  from 
an  arbitrary  trial  wave  function  »p  . 

The  main  purpose  of  this  paper  is  to  show  that  it  is  possible  to  derive 
limits  of  error  for  quantities  of  the  type  !vo(a)  *  r(fO|  >  where  a_  is  an 
arbitrary  point  in  three  dimensional  Euclidian  space,  yQ  is  the  electron  density 
derived  from  the  exact  solution  and  y  is  an  approximate  electron  density  obtained 
from  a  trial  function  ip  .  This  will  be  achieved  in  terms  of  four  quantities;  the 
exact  ground  state  and  first  excited  state  energy,  the  expectation  value  and  the 

mean  square  deviation  of  the  energy.  The  derivation  is  based  largely  on  ideas 

-  3)  — — — — _ 

to  be  found  in  a  paper  by  Kinoshita  .  Kinoshita  estimated  the  error  in  the 

relativistic  corrections  for  a  helium  atom  trial  wave  function.  This  required 

an  estimate  of  the  error  in  the  electron  density  at  the  singularities  of  the  potential. 

It  will  be  shown  here  that  these  results  can  be  generalized  to  apply  on  one  hand 

to  an  arbitrary  point  in  space  and  on  the  other  hand  to  any  atomic  system.  From 

this,  one  will  be  able  to  conclude  that  if  a  wave  function  gives  good  energy  and 

small  mean  square  deviation  for  the  energy,  the  electron  density  will  be  a  good 

approximation  to  the  exact  density,  everywhere,  even  in  those  regions  of  space 

which  from  the  point  of  view  of  energy  may  seem  unimportant. 

It  is  well  known  that  the  state  vectors  which  describe  a  quantum  mechanical 

4) 

system  are  elements  of  the  Hilbert  space  of  quadrstically  integrable  functions  '. 
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These  functions  form  a  Hilbert  space  only  if  the  integration  is  understood  in  the 
sense  of  Lebesque  (von  Neumann  loc.  cit.).  Prom  this  it  follows  that  two  wave 
functions,  which  differ  on  a  set  of  zero  measure  (e.g.  set  of  isolated  points), 
describe  the  same  physical  situation.  Therefore,  it  may  at  first  sight  seem 
surprising  that  one  can  meaningfully  discuss  the  value  of  the  electron  density 
in  a  given  point.  It  will  be  discussed  in  the  next  introductory  chapter  how  the 
ambiguity  on  a  set  of  zero  measure  can  be  removed  if  not  only  i|i  but  also  Hiji 
belongs  to  the  Hilbert  space  of  quadratically  integrable  functions.  This  require¬ 
ment  seems  reasonable  both  from  the  mathematical  and  physical  stand  point. 
The  preparatory  mathematical  considerations  of  the  next  section  are,  however, 
not  strictly  necessary  for  the  understanding  of  the  subsequent  main  discussion. 


2.  PRELIMINARY  DIGRESSIONS 


Let  us  consider  an  atomic  system  with  N  electrons  and  with  nuclear 

charge  Z  .  Let  £.  and  denote  the  position  and  spin  co-ordinates  of  the  i-th 

electron  and  let  x.  stand  for  the  collection  of  r.  and  .  The  Hilbert  space 
1  —l  i 

associated  with  this  system  will  be  denoted  by  L^  (x^,  , .  ,  x^j)  defined  as  the 
set  of  all  functions  ^(x^,  .. .  xN)  which  are  antisymmetric  in  the  variables  xi 
and  satisfy  the  condition 

(<£<£)  4^*iy  )Xv)lZdr.  OO  .  (1) 


Since  the  set  of  Riemann  integrable  functions  do  not  form  a  Hilbert  space  but  only 
a. linear  vector  space,  the  integration  over  £.  is  in  the  sense  of  Lebesque, 
whereas  the  integration  over  the  spin  variable  denotes  summation.  Two 
functions  which  differ  on  a  set  of  zero  measure  (zero  measure  in  3N  dimensional 
Euclidian  space  are  to  be  considered  identical.  The  Hamilton  operator 
H  is  given  in  the  form 

A/ 


H-z-ivt-zf-r-k, 

<-* '  i  =  /  4  i<j  J  ' 

y.  The 


(2) 


where  the  notation  is  self-explanatory, 


operator  H  ,  as  it  stands,  is  not 


yet  well  defined  in  the  mathematical  sense.  Asa  start  one  may  define  H  for 
functions  which  have  second  derivatives  everywhere.  It  was  shown  by  Kato  ^  that, 
provided  the  original  definition  of  H  is  not  too  restrictive,  there  is  one  and 
only  one  self-adjoint  extension.  It  is  necessary  to  make  this  extension  in  order 


-3- 


to  have  a  resolution  of  the  identity  belonging  to  H  i.e.  a  complete  set  of  eigen¬ 
functions.  The  linear  vector  space  on  which  the  self-adjoint  extension  is  defined 

is  called  the  domain  of  H  and  it  will  be  denoted  by  D„  .  D„  is  a  proper 
2  2  H  H 

subset  of  La  i.e.  C  .  As  shown  by  Kato  the  domain  Dpj  consists  of 

all  functions  <j>  t  which  satisfy  the  condition 

r  N  ^ 

J Z-Pi  (5) 

where  4>  is  the  Fourier  transform  of  .  This  may  be  written  somewhat 
loosely  m  the  form 

(^4>ik,P)=  cc?  (4) 

where  K  is  the  kinetic  energy  operator.  It  follows  from  condition  (3)  that  if 

4>  (  D..  then  H<t>  «  L.^  .  that  is,  (H$,  H$)  <  oo  .  Following  the  suggestion  of 
11  7)  A 

Kinoshita  we  require  that  an  acceptable  trial  function  i|i  should  satisfy  con- 

2  2 

dition  (3)  i.e.  we  require  not  only  i]j  c  L^  but  also  Ht|>  €  L^  .  Let  us  now 
define  the  first  order  density  function  y(r^)  as 

He)  =/ <5> 

which  is  i  times  the  probability  of  finding  an  electron  with  arbitrary  spin  at 

the  position  £«.  It  is  a  consequence  of  a  lemma  by  Kato  (loc.  cit. )  that  if  i|i  c  D„ 

H 

y(r,)  is  essentially  continuous  i.e.  it  can  be  made  continuous  by  changing  its 
*  t 

value  on  a  set  of  zero  measure.  Here  we  shall  be  interested  to  estimate  the 
quantity  lro(a)  -  y(a)  |  ,  where  yQ  is  ihe  exact  density  and  a  a  fixed  point. 

The  ambiguity  on  zero  measure  can  now  be  removed  by  the  following  argument. 

In  reality  one  never  measures  the  electron  density  at  a  point  but  only  in  a  small 
volume,  let  us  say  AV  .  It  seems  therefore  reasonable  to  consider  instead  of 
lro(±)  *  r(a)  I  the  quantity 

UvJYoIV^ 'r-Jf  fY(r)dr/' 

Afja  av  ax 

Because  of  the  essential  continuity  of  yQ  and  y  the  result  is  independent  of 
the  way  the  limiting  procedure  is  carried  out  and  of  possible>  discontinuities  on  a 


Incidentally  this  result  is  far  from  trivial.  As  it  is  stated  by  Kato  that 
for  N  >  1  there  exist  functions  which  belong  to  D^,  and  are  still 
essentially  discontinuous,  it  apparently  does  not  hold  for  the  wave  func¬ 
tion  itself. 


sero  measure.  This  is  equivalent  to  redefining  the  densities  so  that  they 
become  continuous  everywhere.  In  the  following  we  shall  thus  be  able  to 
assume  that  the  densities  had  been  chosen  continuous. 


3.  LIMITS  OF  ERROR  FOR  THE  ELECTRON  DENSITY 


Let  ...  x^)  be  an  approximate  atomic  wave  function  which 

satisfies  conditions  (1)  and  (3).  It  can  always  be  written  in  the  form 


V'-0-jtJkyiTVf) 


where  ipQ  is  a  ground  state  eigenfunction  for  the  Hamilton  operator  given  by 
equation  (2)  (the  ground  state  may  be  degenerate),  f  is  a  function  orthogonal  to 
the  ground  state  solutions  and  q  is  a  constant  which  can  be  chosen  real  and 
positive.  If  i]>o  and  ip  are  normalized  then  so  is  f  .  For  typographical 
convenience  we  assume  all  three  functions  to  be  real.  The  constant  q  satisfies 
the  inequality  ^ 


nz<=  ^  -£ 
' 

where  X  =  (Hip,  ip)  and  Eq,  Ej 
state  energies.  The  inequality 


o 


O 


(7) 


are  the  exact  ground  state  and  first  excited 


/o-f)(A%  vt)-ml<vl(AW+HV)h  /)/(8> 


can  easily  be  derived  from  equation  (6)  for  any  linear  operator  A  .  If  A  is 
•elf -adjoint  and  positive  definite  the  identity  A  =  sf&'JK  and  Schwarz's 
inequality  allows  us  to  put  this  in  the  simpler  form 

l(/-if)(%  AV,) V(H>,  AtyjfcAf)  +  (,) 

Let  us  now  put  A  =  6  (a  -  jr ^  ),  where  6(a^  -  r^)  is  the  Dirac  6  function  in 
the  point  a  working  on  the  spacial  co-ordinates  of  electron  i.  One  thus  obtains 

(,0) 

Our  aim  is  to  majoriae  the  right  hand  side  of  inequality  (10)  in  terms  of  Eq, 
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E  ,  X  and  <r  ,  where  r  =  v(H »)/,  Ht(<)  -  X2  .  This  will  be  achieved  in  two 
steps.  First  it  will  be  shown  that  the  unknown  expression  [6(a  -  £j)f.  f]  can 
be  majorized  in  terms  of  (Kf ,  f)  and  (Kf,  Kf)  where  K  is  the 

kinetic  energy  operator.  Secondly,  it  can  be  proved  that  the  integrals  (Kf,  f) 
and  (Kf,  Kf)  on  their  turn  are  bounded  by  an  expression  which  contains  only 
the  afore  mentioned  four  quantities.  Let  us  therefore  consider  [  b(&  -  r^Jf,  f] 
in  detail.  By  the  definition  of  the  6  function 


where  g(r^  4  j  >  £2*  +  a,  tj,  x^,  ...  x^)  and  6(1^)  is  the  Dirac 

6  in  the  origin.  The  identity 

(&&)<},$)-  /&;  <$(r 

/  f  3.  (12) 

=  47 rjr^v,  ■■■,&»))  *T> 


where  dT  a  dr  ^  . . .  di^dtj  . . .  dj^  ,  can  be  proved  by  partial  integration  provided 

j2  2 

f  g  dq)  jdO  jd£2  •••  dr^d^j  ...  d^N  exists  everywhere  except  possibly  on 
£rset  of  isolated  points.  This  leads  to  the  inequality 

r  (13) 

+ j~  (<p*dt  yfdrl Ur  1+  If 


The  first  term  in  the  right  hand  side  expression  may  be  estimated  by  using 
Schwarz's  inequality  as 


On  the  other  hand 


+2  fpr,  fdT 
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which  leads  to 

lffr?VT/* 

and  also 

IfiYYY 

In  similar  fashion  one  may  write 

J~  {fed,  Vt  -  ■  f**')  f lod<  %dT 

since 

i  = v,  !r  rad<  (•  'y*  Y- 


This  gives 


and  therefore 


Substituting  inequalities  (15)  and  (16)  into  (13)  one  obtains 

l 

I  ( kc)  <$  ,y)l*i %{(-^,cj)(9^ ,  7*S )  f 

2 

The  Laplace  operator  V ^  is  invariant  under  translation  by  a  and  therefore 

I  fa-rot,  4)(vxv'4)t 


From  this  it  follows  that 
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•ince  fl  f)  and  C)^ 

^  7j(KA  ^  8u^8t^tut^on  °f  (18)  into  inequality  (10)  gives 


/('-?)  Yj*> -^3; /  <2^{^  fkWyjll KM] 

*}\ fwT)  TjUKt! . 


/ 


-r- 


-f-1  u 
7 r  v 


(19) 


The  first  part  of  our  programme  has  now  been  accomplished,  si  (Kf,  i) 
and  || Kf  ||  are  the  only  unknown  quantities  in  the  right  hand  side  expression. 
These  can  be  majorized  as 


yllmi  < ]f^y--(l-rj2-) £f  +  \fc rj lA W, V 


) 


(20) 

(21) 


where  a  is  a  numerical  constant  given  by 


0L~J.(4Z2--tlZ.(2.-l)(U-i)  f  N(N-lf) .  (22) 

For  details  the  reader  is  referred  to  Appendices  I  and  II.  If  one  introduces 
the  notation 

C-f%(z  l/?7  *-£°)  +tI‘£°  )  ■ 

y  '(23) 


it  will  follow  from  inequalities  (19).  (20)  and  (21)  that 

I  0  -c  c  <24> 

which  leads  to  the  final  formula 

/  Y/ai-fi a) I  <  j—z  (i  ifc  l/itz)  •+  C  +  *fY<&)) 

where  C  is  given  by  equation  (23)  and  tj  satisfies  inequality  (7). 


(25) 
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We  may  thus  conclude  that  the  error  in  the  electron  density  !iro(a)  * 

•  y(a)  |  ,  at  any  given  point  a  ,  is  bounded  by  an  expression  which  depends 
besides  r(a)  only  on  the  quantities  X,  or,  Eq  and  q  .  Inequality  (7)  gives 
an  upper  bound  for  q  whereas  X  and  < r  can  be  calculated  from  the  wave  func¬ 
tion  tj<  .  As  for  Eq  and  Ej  one  may  either  use  experimental  values  or  else 
it  is  possible  to  give  upper  and  lower  bounds  by  known  methods. 

It  also  follows  from  formula  (25)  that  if  a  sequence  of  trial  functions 
converges  to  the  exact  solution  i.e.  ||  -  «|»  ||  -♦  0  and  if  in  addition 

||  (H  -  X  ||  -*  0  ,  the  convergence  to  the  electron  density  is  uniform. 

It  is  to  be  noted  that  these  arguments  can  easily  be  generalized  to  apply 

to  a  molecular  system.  The  only  difference  in  the  final  formula  would  be  that 
n 

now  Z  =  £  Z  ,  where  n  is  the  number  of  atoms  in  the  system. 
a=  1  a 

A  final  remark  concerning  the  spin  density  should  perhaps  be  made- 
The  spin  density  at  a  point  a  may  be  defined  as 

~J &(T'-£)'ty>(X'n'~>xJS2i  ^{Kiy..)^u)dbcv..cbcN 


where  the  spin  matrix  works  only  on  the  co-ordinates  of  electron  1.  Si 

II  sz  II  =  i  it  can  be  shown  by  using  Schwarz's  inequality  that 


/ £  (oZ  fifte)  [c rfyiOz) ) 

where  y(a)  is  again  the  total  electron  density.  The  spin  density  plays  an 
important  role  in  the  di8CUB®i°n  of  the  hyperfine  interaction  and  the  above 
formula  might  be  of  use  for  testing  the  reliability  of  very  accurate  calculations. 


4.  LIMITS  OF  ERROR  FOR  ATOMIC  FORM  FACTORS 

The  form  factor  is  defined,  apart  from  an  unimportant  numerical  factor, 
as  the  Fourier  transform  of  the  electron  density  Let  therefore  ro(k)  and 
r(k)  denote  the  Fourier  transforms  of  T3(£)  and  v(r)  •  By  substituting 
A  =  — -*77T  e  *  —  ‘  —  in  inequality  (8)  one  obtains 

iW)Y&-fa>h  ^(is£*'cw^ 


(26) 
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and  because  of  Schwarz's  inequality  this  leads  to 

l(i-yl)  YJJs) -  Wi &  vaf)*/* 

since  tj*  and  f  are  normalized.  This  gives 

+77a/fffei/)  <27> 

Moreover,  since  |y(k)  |  <  —  -■ -T  / y(r)dr  =  -*r-  one  also  has  the  formula 

(2ir  l'**  ~  ~  (2it)7a 

/K& -/(&!<£ . 

The  last  two  inequalities  might  be  useful  for  estimating  the  accuracy  of  the 
form  factors  used  in  crystallography. 


5.  NUMERICAL  EXAMPLE 

Here  we  shall  illustrate  the  use  of  formula  (25)  by  applying  it  to  the 
helium  atom.  Since  for  helium  Z  =  2  and  N  =  2  we  have 

Y&>)  <28> 

where 

( 1 / 1/(40^+  Q-ej )\ 

Here  we  made  use  of  the  fact  that  r\“  is  bounded; 

_  2 

Table  I  lists  the  values  of  C\  2 s/C7  and  ti  for  a  set  of  wave  func- 

1  max 

tions  containing  6,  18  and  38  variational  parameters.  For  the  exact  ground  state 
and  first  excited  energies  Eq  =  -2. 903  725  a.u.  and  Ej  «  -2.  146  a.u.  were 
used  (reference  3). 

It  seems  that  inequality  (2  8)  gives  sensible  error  limits  at  least  for  very 
accurate  wave  functions. 
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2 

TABLE  I.  The  calculated  values  of  2n/Ct,  C'  and  ri  for  helium  atom 

max 

wave  functions. 


No.  of 

parameters 

-X 

2 

<r 

2 

n 

1  max 

2*/C1 

C' 

6 

2.903  24 

0.016  90 

0.000  485 

0.40 

0.040 

18 

2.903  715 

0.000  922 

0.000  013 

0.0602 

0.000  906 

38 

2.903  722 

0.000  115 

0.000  004 

0.032 

0.000  246 

The  value  8  of  -X  and 


2 

<r 


were  taken  from  Table  II  in  reference  7. 
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APPENDIX  I 


Upper  bound  for  f) 


Here  we  shall  sketch  the  proof  of  inequality  (20)  in  section  3.  From 
equation  (6) 


it  follows  that 


=  3-  E0+ifE< 


which  may  be  written  as 


where 


V-  -1 4  +Z-  r", 

-  Uj  ! J 


.  Clearly 


and 


(vw>-z(iu4j 

also  iml 


(29) 


=■  %  -  o )  (3°) 


im  /  1 

For  details  Bee  Kinoshita  (1959)  loc.  cit.  The  substitution  into  equation  (30) 
leads  to 

rfOctf  )-rf7.iw  1 /(ict.t) «  ^  +  V*E° 


which  can  be  solved  for  tj n/ (Kf,  f)  to  give  the  desired  formula 


v  <  z  V¥ i- 
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APPENDIX  II 
Upper  bound  for  q  ||  Kf  || 


As  K  =  H  -  V  . 


r}HW<vMI-fHVU 

while  from  equation  (6)  it  follows  that 

l)tfK/z{ra+V-(i-V')E;]* 

where  <r  =  ||  (H  -  X  ) 4*  II  •  This  gives 

7  IK#  «  {<r2+V-(/-  i]z)E!  } 7  //!/■///.  (31) 

It  remains  to  show  that  r)  ||  Vf  ||  is  bounded 

( vi  m  =  ( v%  4. ) = a  -zz*c )%  i ) = 

=  2  A + 

-  ■* Z  A/fA/-?j ( j  £ {tl)-z  A K v-0( A/-2) ^  {) 

since  f  is  antisymmetric  in  the  electron  co-ordinates.  By  dropping  the  term 
-ZN(N  -  l)(N  -  2)(~  -pi—  f,  f)  which  is  always  negative  and  applying  Schwarz's 
inequality  to  the  thre^  an^  four  body  terms,  one  obtains 

(vm) 

+zzN(N-i)(t'  iw--zz  »(* /-o(i  tw. 


The  last  two  terms  cancel  partly.  This  is  a  consequence  of  the  inequality 
rl2^  rl  +  r2  *rom  w^tch  it  follows  that  ■—  -  —  — —  -  ~  t~—  0 


1 


'l  r  12 


12 


-13- 


and  therefore 


(yi,  1/-AK2,  A l(r,'t)^,£)  +  rui) 

^  Z(z.-i)N(V-,)  (^l  ±f) 


+ 


It  was  shown  by  Kino shita,  loc.  cit.  (1959).  that 

r,{)  s  4  (-^44,4)  , 


S=2(~VXJ) 


and  so  one  obtains 

(v{,W)^cc(K(,f)  HV4H 

wh=r.  ct--JL(4Z,^2z.(z.-l)(h/-i)  -f-  N(N-lf) 

Substituting  this  into  inequality  ( 3 1 )  the  desired  formula 

7 / <{<T+ Jf-O-YlE?}  A  i/S  7  \/(K4,4) 

is  obtained  as  it  was  shown  in  Appendix  I  that  q«J"(kf,~f)  is  bounded.  This  is 
a  generalization  of  Kinoshita's  formula  for  the  many-electron  case. 
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